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ABSTRACT 

We present one-dimensional models of the hot gas in dark-matter halos, which both 
predict the existence of cool cores and explain their structure. Our models are directly 
applicable to semi-analytic models fSAMs) of galaxy formation. We have previously 
argued that filaments of cold (~ 10 K) gas condense out of the intraclustcr medium 
(ICM) in hydrostatic and thermal equilibrium when the ratio of the thermal instability 
timescale to the free-fall time tn/iff falls below 5-10. This criterion corresponds to an 
upper limit on the density of the ICM and motivates a model in which a density core 
forms wherever iTiAff < 10- Consistent with observations and numerical simulations, 
this model predicts larger and more tenuous cores for lower-mass halos — while the core 
density in a cluster may be as large as ~ 0.1cm" 3 , the core density in the Galactic 
halo should not exceed ~ 10 _4 cm~ 3 . We can also explain the large densities in smaller 
mass halos (galactic 'coronae') if we include the contribution of the central galaxy to 
the gravitational potential. Our models produce a favorable match to the observational 
X-ray luminosity-temperature (Lx — Tx) relation. For halo masses <, 10 13 M Q the core 
size approaches the virial radius. Thus, most of the baryons in such halos cannot be in 
the hot ICM, but either in the form of stars or in the form of hot gas beyond the virial 
radius. Because of the smaller mass in the ICM and much larger mass available for star 
formation, the majority of the baryons in low mass halos (<; 10 13 M Q ) can be expelled 
beyond the virial radius due to supernova feedback. This can account for the baryons 
'missing' from low mass halos, such as the Galactic halo. 

Key words: galaxies: clusters: intracluster medium; galaxies: halos. 



1 INTRODUCTION 

In the currently favored ACDM cosmology, massive halos 
form due to mergers of smaller dark matter halos. The over- 
all structure of the halos is simply an outcome of gravita- 
tional inter actions between dissip ationless dark matter par- 
ticles (e.g., ISpringel et al. 1 120051 ), and once the dark mat- 
ter halo relaxes the density tak es a simple self-similar form 
(|Navarro. Frenk. fc White 1 19971 ). While the evolution of dark 
matter halos is simple, at least in principle, the evolution of 
baryons residing in these halos is much more complicated. In 
addition to gravity, the baryons are affected by complex, dif- 
ficult to model processes such as radiative cooling, feedback 
heating, and turbulence. 

Observationally it is well known that the distribution 
of baryons in the form of hot gas and condensed phases 



(stars, molecular and atomic gas) is a sensitive function of 
halo mass. Baryons in the most massive halos (clusters and 
groups) are mainly in the form of the hot, X-ray emitting 
plasma. The baryon fraction in the hot phase decreases with 
a decreasing halo mass, but the fraction in stars increases 
to keep the tot al baryon fraction roughly constant for ha - 



los £ 10 13 M (|Gonzalez et al. ll2007l ; TGiodini et al. 1120091 ) 
While massive halos seem baryonically 'closed', smaller ha- 
los, inclu ding our own M i lky Way, have lost mo st of their 
baryons l|Dai et al. I l201Ct IMcGaugh et al. I |201Ch . Observa- 
tionally, the halo mass seems to be the primary variable deter- 
mining the halo baryon and ste llar fractions, irrespective of 
AGN activity or star fo rmation (|Anderson fc Bregman Ifeoid ; 
IMcGaugh et al. Il20ld ). This remarkable correlation between 
the halo mass and stellar/baryonic mass provides an impor- 
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tant clue to understanding galaxy formation. A related ob- 
servational result is that the smaller halos have lower X-ray 
luminosities (Lx) than self-similar models which assume a 
similar gas density /temperature profile for all halos. In this 
paper we present models which explain these observations in 
terms of the interplay of cooling, feedback heating, and ther- 
mal instability in halos. 

Cooling and feedback play a fundamental role in galaxy 
formation. Both cooling and star formation are efficient 
if the cooling time at the virial radius is shorter than 
the dynamical/ free- fall time (e.g., iRees fc Ostriker I 1 19771 : 
IWhite fc Rees I fl978). Even when the cooling time is longer 
than the free-fall time, it can be shorter than the halo age 
(~ Hubble time). Thus, in absence of additional heating, the 
inner radii of halos affecte d by cooling a re expected to form 
massive cooling flows (e.g.. IFabian 1 1994). However, both the 
absence of soft X-ray lines (e.g., Peterson et al. 1120031) and th e 
lack of significant star formation (e.g.. lO'Dea et al. I|200S| ). 
indicate that cooling in cluster cores is largely suppressed. 
There is growing evidence that a ctive galactic nuclei (AGN) 
jets, bubbles, and cavities (e.g.. iBirzan et al. 1120041 ) blown 
by the central supermassive black holes (SMBHs) in clus- 
ter cores play a crucial role in suppressing cooling flows (see 
iMcNamara fc Nulsen |2007 | for a review). S upernova feedback 
in lower mass halos (e.g- lEfstathiou I feoOO ) appears sufficient 
to prevent runaway cooling of the hot gas. Moreover, feedback 
heating is globally stable because t he condensation of cold ga s 
leads to enhanced heating (e.g., see lPizzolato fc Soker 1 ,2005). 
Understanding the hot gas thermodynamics is essential to ex- 
plain the non-self-similarity of baryons in halos, as we argue 
below. 

Motivated by observations which suggest balance 
between heating an d coolin g, in a series of papers 
dMcCourt et al. |2012| . hereafter lPaper itlSharma et al. l2012l 
hereafter iPaper HI ) we studied the role of local thermal insta- 
bility in governing the structure of the pressure-supported 
hot gas (~ 10 6 — 10 8 K) for virialized halos in global thermal 
balance. Our idealized simulations, which impose a balance 
between radiative cooling and heating, predict that cold fil- 
aments condense out of the hot phase if, and only if, the 
thermal instability timescale (£ti) is shorter than approxi- 
mately ten times the free-fall time (tff)Q This criterion gives 
an upper limit on the hot gas density in cores of different 
halo masses. For gas densities larger than this, cold filaments 
condense out of the hot phase. The cold overdense filaments 
are heavier than the surrounding hot gas, and thus fall to- 
ward the center and power AGN/supernovae. Thus, even if 
feedback heating can balance cooling in the cluster core, the 
physics of local thermal instability gives a rough upper limit 
on the density of the hot gas. The end result is that the core 



1 Cold gas of any kind is referred to as filaments in this paper; 
numerical simulations and observations suggest that the cold gas 
is filamentary. For a heating rate per unit volume which is indepen- 
dent of density, the thermal instability timescale is roughly equal 
to the cooling ti me; for definitions of various timescales see Eqs. 
8-10 in lPaper ij . 



density is reduced until the criterion tTl/ig > 10 is satisfied 
for the hot gas. 

Our cores with tTl/tg ~ 10 agree with the ob- 
servations of cool core clusters. Our results are robust 
as long as the timescale for secular cooling/heating of 
the ICM is longer than the radiative cooling time; i.e., 
if the halo i s in ro ug h therm al balance. The reader 
should consult IPaper 1 fc IPaper il for more details. Recent, 
more realistic cluster simulations with AGN jet feedback 
(|Gaspari. Ruszkowski. fc Sharmal [2012) roughly agree with 
the findings of our idealized simulations; in particular, ex- 
tended multiphase gas condenses out of the ICM and leads to 
enhanced feedback when our r-Ti/iff < 10 criterion is satisfied. 
However, as expected, there are differences in detail; e.g., the 
presence of angular momentum supported cold torus at small 
radii. 

Based on our in/is criterion for the formation of multi- 
phase gas, it is straightforward to build one-dimensional hy- 
drostatic models for the stru cture of the hot gas in virialized 
halos (see IVoit et al. 1 12002| for a similar approach). We as- 
sume a universal gas density profile at large radii unaffected 
by cooling/feedback, and we build hydrostatic profiles inward 
until the radius where tTi/te is smaller than the threshold 
for forming cold filaments. An isentropic core is introduced 
within that radius, a s suggested by our nume rical simulations 
and by observations l|Cavagnolo et al~ll20 09) . We explore our 
models by changing several parameters, the outer gas density 
profile, tTi/ta threshold, a constant <Ti/tfl core, inclusion 
of stellar gravitational potential, etc. Our results show that 
lower mass halos have larger but lower density cores because 
of short er cooling t imes. This is consistent with group obser- 
vati ons JSun 1120121 ) and with numerical simulations (see Fig. 
9 in IPaper III) . The X-ray luminosity-temperature (Lx — Tx) 
relation derived from our models is also consistent with obser- 
vations. Moreover, we do not expect much redshift evolution 
of the hot gas in virialized halos because the tTi/tg criterion, 
and hence the core density, is essentially independent of the 
redshift. 

We estimate the feedback efficiency required to balance 
radiative cooling (Lx) in halos with different masses. We find 
that large feedback efficiencies (0.001 — 0.01) are required for 
massive clusters. However, for lower halo masses the core is 
extended and has a lower density. Assuming that most of the 
baryons that are absent from the ICM are channeled into star 
formation and accretion, implies that a smaller feedback ef- 
ficiency is required for lower mass halos such as groups. For 
halos of even smaller mass (<; 10 13 Mq) the required feedback 
efficiency is extremely small (< 10~ 7 ), 10— 100 times smaller 
than what is available due to supernova feedback. This means 
that most of the baryons in such halos will be overheated 
and will be pushed beyond the virial radius. Therefore, for 
lower mass halos, such as our own Galaxy, most baryons 
which should have been accreted in the absence of feedback, 
are 'missing'. These 'missing' baryons are in the extended 
corona which is prone to forming multiphase gas required 
to ex plain quasar absor ption lines observed close to galaxies 



to explain quasar absorp 
fe.g Jchen et al. Il2010l ). 
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This paper is organized as follows. We describe our one- 
dimensional models in Section 2. In Section 3 we derive the X- 
ray luminosity-temperature (Lx—Tx) relation from our mod- 
els and compare it with observations. Section 3 also presents 
the implications of our models for the 'missing' baryons prob- 
lem. We conclude in Section 4 with astrophysical implications 
of our models. 



ONE-DIMENSIONAL MODELS FOR HOT 
GAS IN HALOS 



We c losely foll ow th e compr ehensive paper by Voit et al. I 
|2002| (hereafter |V02| ; see also iFukugita fc Peebles 1 120061 ) to 
build our one-dimensional models. Nevertheless, there are 
significant differences. The key difference is that they apply 
an entropy floor within the radius where the cooling time is 
shorter than the cluster age, but we apply an entropy floor 
(we also try a floor on tTi/te, as described later) where our 
tTl/tg < 10 threshold is satisfied. The physical p ictur es cor- 
responding to the two criteria are different. The V02 model 
posits that the X-ray gas with a cooling time shorter than 
the cluster age either cools to the cold phase or is heated by 
feedback such that the cooling time is longer than the cluster 
age. Theoretically, however, steady-state cooling flows main- 
tain cooling times much shorter than the Hubble time. The 
cooling scenario is further disfavored as it suffers from the 
overcooling problem. The core entropy obtained by equating 
the cooling time to the cluster age (~ 5 Gyr; e.g., Fig. 1 in 
IVoit fc Bryan ifaoOll ) is much larger than wha t is observed in 
cool core clusters fe.g.. lCavagnolo et al. l[2009h . and there are 
also non-cool core clusters with cooling times longer than the 
cluster age. We note, of course, that the lV02i model seems to 
satisfactorily describe the non-cool core cluster density pro- 
files and the resul ting L x — Tx relation. 

In contrast to|V02J, our threshold is based on the physics 
of local thermal instability in cluster halos. The central as- 
sumption of our models is that heating and cooling roughly 
balance globally and regulate the core to an approximate ther- 
mal equilibrium. While we do not understand the details of 
how the ICM is heated, the assumption of rough thermal bal- 
ance is robust and observationally motivated for cool core 
clusters. The lack of massive cooling flows suggests that the 
heating rate should be at least as large as the radiative cool- 
ing rate. Likewise, heating cannot be too dominant; other- 
wise the cold gas would be totally e vaporated (see Fig. 13 in 
Sh arma. Parrish. fc Quataert 1120101 ). Thus, cluster cores are 
expected to be in rough thermal balance over a few cooling 
times. 

The physics of local thermal instability in the ICM (we 
refer to the hot gas in all virialized halos, even groups and 
individual galaxies, as the ICM) in thermal balance intro- 
duces a new parameter tTi/iff which is absent in earlier one- 
dimensional models. The thermal instability occurs in two 
regimes depending on txi/iff- If this ratio is smaller than a 
critical threshold (~ 10) then local thermal instability leads 
to multiphase gas; an overdense blob cools to the stable tem- 



Thus, no multip hase gas is expected if tTi/iff J> 10. We ar- 
gue in lPaper III that, since extended multiphase gas leads to 



perature and falls toward the cluster center on ~ free-fall 
timescale. For a larger in /iff, the overdense blob responds 
to gravity as it is cooling. Shear generated between the in- 
falling overdense blob and the background hot gas is able 
to mix it before it can cool to the thermally stable phase. 
no multip 
1 [Paper III 1 

very high accretion rates, cooling and feedback in cool cores 
should self- regulate to the critical threshold where in ~ lOtft . 
Our in /iff criterion results in a realistic core entropy for cool 
core clusters and performs well as a predic tor of extended 
cold gas in cluster cores (e.g., see Fig. 12 in lMcCourt et al."1 
2012). The physical basis of our models and good quantitative 
match with observations give us confidence in the predictions 
of our one-dimensional models. 

We point out that, though our criterion on the ratio 
i n /iff i s very different from the criterion on the cooling time 
in |V02| for group and cluster masses, the two models yield 
similar results for Galactic mass halos. This is because we 
find that the core size in the Galactic halo is ~ the virial 
radius, where the dynamical time is ~ l/10 th of the Hub- 
ble time. Thus, for 10 12 solar mass halos, our criterion that 
in /iff ~ 10 is equivalent to the as sumption that the c ooling 
time roughly equals the halo age. iMaller fc Bullockl l|2004h 
explored models of the Galactic halo in which the cooling 
time at the virial radius is matched to the halo age; our re- 
sults in this mass range agree closely with theirs. The long 
cooling time at the virial radius also signals the breakdown 
of our model, however: if the cooling time is comparable to 
the age of the system, our assumption of thermal equilib- 
rium need not hold. Thus, our results for halo masses below 
~ 6 x 10 12 Mq should be viewed more as an extrapolation 
from higher masses than as solid quantitative predictions. 



2.1 Recipe for Constructing One-dimensional 
Models 

We construct one dimensional gas profiles such that the 
ICM is in hydrostatic equilibrium and everywhere satisfies 
our criterion in/iff J> 10. In determining hydrostatic equi- 
librium we neglect gravity due to gas and stars for most 
models, and assume a fixed N FW dark matter potential 
(|Navarro. Frenk. fc White Hl997T l. with the concentration pa- 
rameter (c = r2oo/r s , where r2oo is the virial radius within 
which the average density is 200 times the critical density 
and r s is the scale radius) adjusted according to the halo 
mass. The relation between the concentr ation parameter (c) 
and M200 (mass within r2oo) is based on |vo2 and references 
therein (see Table [1]). Our results are only weakly sensitive to 



Following we construct an 'unmodified' profile in 

which the gas density is specified and other thermodynamic 
variables are calculated using hydrostatic equilibrium. The 
'unmodified' density profile agrees with the observed ICM 
density at large radii and has a baryon fraction equal to the 
cosmic value. Thus, the total gas mass in the 'unmodified' 
halo equals fbM2oo, where f b = 0.17 is the universal baryon 
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fraction. Although observations suggest that smaller halos are 
baryon poor, for simplicity we use a fixed baryon fraction for 
all 'unmodified' halos. As we show in Section 13.21 the final 
baryon fraction is something we can calculate using our mod- 
els. Our results are only weakly dependent on the assumed 
initial value of ft. 

The effect of cooling and heating on the ICM is in- 
corporated by altering the initial profile so that it satisfies 
in /iff > 10 everywhere. This process amounts to introduc- 
ing a low density core in the center of the cluster. Since this 
prescription is not unique, we try two different core profiles: 
an isentropic (T/p 2 ^ 3 = constant) core; and a 'core' with a 
fixed tTi/tg. Each constraint, along with hydrostatic equilib- 
rium, then uniquely specifies the structure of the core. In each 
case, the atmosphere outside the core is unmodified. 

Observations suggest that an isentropic core is a good 
approximation for the inner radii o f the ICM, for both coo l 
core and non-cool core clusters (see ICavagnolo et al. 112009^ 1 . 
A 'core' with iTi/iff equal to the critical value for forming 
cold filaments represents an upper limit on the density of the 
hot phase at all radii within the core. If the density exceeds 
this limit, a large mass of cold gas will condense out of th e 
hot phase, leaving behind a core with iTi/iff ~ 10 (jPaper I J ). 

Since the gas density profile, even in cluster outskirts, is 
not universal (due to mergers, accretion hi story, cosmologi- 
cal filaments, etc.; e.g.. ICroston et al. 112008^ 1 . we try different 
'unmodified' gas density profiles of a generalized NFW form 

= Nc m 

9 r/r.(l + r/r,)(— D' 1 ' 

where jV c is a normalization constant such that 



4ivr pdr = fbM% 



and s is the asymptotic (r — > oo) gas density slope dlnp/dlnr 
[s = —3 for an NFW density profile). 

The last parameter needed to specify our models is a 
boundary condition on one of our thermodynamic variables. 
We choose the pressure at V200 as the outer boundary condi- 
tion. The outer pressure is specified as P200 = fc s (n)fc.B T200, 
where k s is a constant, (n) = 3/i,M2oo/47r^m p r2oo (m is 
the mean molecular weight, m v is the proton mass) is the 
average gas number density within r2oo, and fcflT2oo = 
GM2ooA tm p/2?"20o is the virial temperature at V2oo- The con- 
stant k s is smaller than unity since the number density at 
r"2oo is smaller than the average density (n); moreover, k s 
should be smaller for a steeper initial gas density profile. We 
choose k a = 0.6, 0.32, 0.15 for s = -1.75, - 2.25, - 3 
respectively. These values ensure that a power-law entropy 
profile is obtained at large radii for 'unmodified' profiles in 
hydrostatic equilibrium. The results are quite insensitive to 
the exact value of k s . 

We calculate our entropy/temperature profiles starting 
at r2oo and moving in by imposing hydrostatic equilibrium 
and using the ICM density given by E q. [T] We calculate the 
iTi/iff ratio (see Eqs. 8 -10 in Paper II for defini tions) at each 
radius. We use t he fit inlTozzi fc Normaiil J200ll) (Eq. B4, Ta- 
ble 3) based on ISutherland fc Dopita I l|l993h for the cooling 



function corr esponding to one-third of the solar metallicity 
(see Fig. 1 in lSharma. Parrish. fc Quataert Il2010h . If in /iff 
at some radius becomes smaller than the critical value (chosen 
to be 10 for most of our models) we introduce a 'core' within 
that radius. We use the same critical value of tn/iff irrespec- 
tive of the halo mass. This assumption is justified if halos 
are in rough thermal balance and the microphysical heating 
mechanism is similar for all of them. 

For smaller halos (< 10 13 Mq), ixi at the virial radius can 
be shorter than lOtg because of efficient line cooling. In these 
cases, we adjust the density and the temperature at the outer 
boundary such that tTi/iff there equals the critical value and 
the pressure equals the imposed value at the outer bound- 
ary. Since the results of our one-dimensional models in such 
cases depend on the outer boundary conditions (e.g., whether 
we keep the outer pressure/temperature/density fixed), the 
structure of the hot plasma in lower mass halos can only be 
studied in detail by cosmological numerical simulations. How- 
ever, most of our results stem from the trend of core size with 
halo mass, which we expect to hold even for the lower mass 
halos. It is important to note that the limit on iTi/iff (~ 10) 
is a lower limit which should apply to cool core halos; the 
minimum tTi/iff for gas can be larger than 5 — 10 in non- 
cool core systems because of overheating by feedback and/or 
major mergers. 

Here we briefly summarize the various steps involved in 
constructing our one-dimensional models: 

(i) Assume an 'unmodified' gas density profile of a gener- 
alized NFW form (Eq. [TJ|, such that the total gas mass is ft 
times the halo mass. 

(ii) Specify the pressure at V2oo as an outer boundary con- 
dition and construct 'unmodified' profiles by assuming the gas 
density profile given by Eq. [T]and imposing hydrostatic equi- 
librium in an NFW potential ( we also include gravity due to 
the central galaxy in some models); choose P200 such that the 
entropy profile close to the virial radius is a power-law. 

(iii) Calculate iTi/iff at each radius moving inside from the 
outer boundary; introduce a 'core' with a constant entropy 
or a constant in /is within the radius where iTi/iff becomes 
smaller than the critical value for cold gas condensation. 

(iv) Impose hydrostatic equilibrium within the 'core' by 
assuming a constant entropy (or ii-i/iff) and obtain the 'mod- 
ified' profile. 



3 RESULTS 

Table[T]contains our models for halos ranging from the Galac- 
tic mass (10 12 Mq) to most massive clusters (3 x 10 15 Mq). 
This table lists various quantities calculated from our con- 
stant entropy and constant tTi/iff models with the critical 
iTi/iff = 10. The table lists the central electron density, the 
central entropy, the X-ray luminosity-weighted temperature 
(Tx), and the X-ray luminosity (Lx)- Table Q] also lists the 
mass dropout factor fd, the fraction of the baryonic mass 
which is removed from the hot ICM when we 'modify' the 
core. This mass can either go into fueling star formation and 
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Figure 1. The electron number density (n e ) as a function of radius 
(scaled by r2oo) f° r different halo masses; the asymptotic gas den- 
sity slope s = —2.25 for these models. The profiles for an entropy 
core with txi/iff ^ 10 are shown by the solid lines and the profiles 
for a fixed tTlAff = 10 'core' are shown by the dotted lines. The 
halo masses, in the order of decreasing core density, are 3 X 10 14 , 
6 x 10 13 , 10 13 , 6 x 10 12 , 3 x 10 12 , and 10 12 M Q . 



black hole growth, or can be pushed outside T2oo because 
of strong feedback heating (in this sense the term 'dropout 
factor' is misleading because it includes hot gas beyond the 
virial radius). In addition, the table lists an estimate of the 
mass dropout rate Md, the ratio of the mass absent from 
the hot ICM and the typical halo age £ age . The halo age is 
assumed to be 5 Gyr for all halos; merger rate, and hence 
ta, ge , is very weakly dependent on the halo mass; e.g., see 
iFakhouri. Ma. fc Bovlan-KolchirTI l201fj| . From M d , we esti- 
mate the required feedback heating as the ratio of the X-ray 
luminosity and MdC 2 , e rC q = Lx/MdC 2 . 

Figure [T] shows the electron number density profiles 
for different halo masses using our constant entropy (solid 
lines) and constant tTi/ts (dotted lines) core models; the 
critical tTi/ta is 10 and the asymptotic gas density slope 
s = dlnp/dlnr = —2.25. The halo mass ranges from that 
of the Galactic halo (10 12 M ) to a cluster (3 x 10 14 M©). 
Figure [T] shows several trends with the halo mass. Most im- 
portantly, the core electron density is smaller and the core 
size is larger for smaller mass halos. Cool core clusters with 
mass > 10 14 Mq have similar scaled density profiles according 
to our models, especially outside 0.05r2oo- This is confirmed 
by observations of relaxed clusters with temperatures > 3.5 
keV (see Fig. 13 in lMaughan et al. II20I2I ); density in smaller 
cool core halos is lower, in agreement with our models. The 
constant tTi/tg cores give a much larger central density com- 
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Figure 2. Electron number density as a function of radius for 
cool core ACCEPT groups and clusters (core entropy < 30 keV 
cm 2 ; ICavagnolo et aL~ll2009l) , and our models including the BCG 
potential. ACCEPT data is averaged in 5 temperature bins: < 
1 keV (< 3.5 X 10 13 M ; grey filled circles); 1-2 keV (3.5 X 
10 13 - 1.1 X 10 14 M Q ; red stars); 2-4 keV (1.1 - 3.6 X 10 14 ; blue 
triangles); 4 — 6 keV (3.6 — 7 X 10 14 Mq; green squares); 6 — 8 
keV (7 — 11 X 10 14 Mq; violet hexagons). The temperature range 
is converted into mass range using the best fit M200 — Tx relation 
for our core entropy models in Table [T] M200 ~ 3 X 10 14 Mq 
(Tx/3.6 keV) 5 / 3 . The theoretical profiles are shown for halo masses 
oflO 13 , 6xl0 13 , 3xl0 14 , 6xl0 14 , 10 15 M Q , in order of increasing 
density. 



pared to the constant entropy cores, especially for smaller ha- 
los, because tg decreases toward the center and the density 
has to increase to keep frri/iff constant. 

It has been observed that low mass halos (groups and 
elliptical galaxies) host cool 'coronae' ins tead of co ol cores 
which are observed in massive clusters (e.g JSun 1 20091 ) . These 
coronae have been likened to mini cooling flows which feed the 
central SMBHs and power radio jets. Instead of flat density 
cores, these coronae have density profiles steeply rising toward 
the center. We can explain the coronae in the context of our 
models if we include the gravitational potential of the cen- 
tral brightest cluster galaxy (BCG), in addition to the NFW 
potential. We model the gravitational acceleration due to the 
BCG as given by (in cgs units) 



3BCG 



(r/lkpc) 



3.206 x 10- 7 



+ 



(r/lkpc) 



1.861 x 10" 6 



-l/s 



(2) 



with s = 0.9, and as show n in Figure 1 of 
iMathews. Faltenbacher. fc Brighenti I (2006). We use the 
same BCG potential for all our halos. While this assumption 
is not correct in detail, it captures the essential feature that 
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gravity at the smallest radii is dominated by the BCG and 
not by the NFW potential. 

Figure [5] shows a comparison of our models (which in- 
clude the BCG potential) with observations. The observations 
agree with our models in that the smaller halos have lower 
densities compared to self-similar modles. While we make the 
simplifying assumption that the gas density profile outside 
the 'core' is the same for all halos, observations show that 
smaller halos have shallower density profiles. This is because 
realistic feedback can affect the gas even beyond the core, es- 
pecially for the lower mass halos where the required feedback 
efficiency to balance cooling is quite small (see section |3T2"1 and 
Fig. [4]). A comparison of Figures [T] and [2] shows that the core 
density for iTi/tfl = 10 models is much larger when the BCG 
potential is included. The effect is much more pronounced for 
smaller halos in which the BCG gravity dominates the NFW 
gravity even at relatively large radii. Thus, by including the 
BCG potential we are able to reproduce the observed steep 
rise in density toward the center in our low mass halo models. 

Most of our results stem from the fact that the smaller 
halos have lower core densities; in Section 13.11 we study its 
implications for the Lx — Tx relation, and in Section \3. 21 we 
apply it to the missing baryons problem. 

3.1 L x - T x Relation 

Figure [T] illustrates that the physics of thermal instability in 
halos results in big ger bu t lower density cores for smaller ha- 
los (cf. Fig. 9 in lPaper ij) . This reduces the X-ray luminosity 
of low mass halos below the predictions of gravitational self- 
similar models. Moreover, Table[T]shows that the core entropy 
is nearly independent of the halo mass above > 3 x fO 13 M Q , 
also in contrast with self-similar models. Thus, our models 
explain the observed s t eepening of the Lx — Tx relation 
(|Evrard fc Henry] Il99ll ; iBrvan Il2000l and references therein ) 
and excess entropy (e.g.. |Pom~nan. Cannon, fe Navarro II 19991 ) 
in lower mass halos. 

We can understand this result intuitively as follows. For 
self-similar evolution of hot gas in halos, we expect the gas 
density to be the same for all halo masses. This assumption 
leads to L x tx n 2 A{T)rl a0 oc T 3/2 A(T), which goes like T 2 for 
free- free cooling (A oc T 1 ' 2 ). On the other hand, holding <ti <x 
T/nA(T) constant for all halos (assuming a similar tg for dif- 
ferent halos) gives n oc T/A(T), which is oc T 1//2 or M^qq 
for free-free cooling. Thus, our models predict lower densities 
for smaller halos, and hence Lx oc n 2 A(r)rfoo oc T 3 for mas- 
sive halos, in line with observations. Above arguments assume 
that most of the X-ray luminosity is contributed by the dense 
core and that the scaled core size is the same for all halos; 
as Figure [1] shows, these assumptions are not satisfied quan- 
titatively. The extra luminosity contributed by a bigger core 
for smaller halos is partly compensated by a smaller density 
because tg (roughly oc r 1//2 ) at larger radii is longer. Later 
we quantitatively discuss the Lx — Tx relation obtained from 
our models. 

At halo temperatures below 1 keV, where line cooling 
becomes important and the cooling function starts to increase 



with a decreasing temperature, our ansatz of keeping a fixed 
tTi/ta yields an even smaller core density and an even steeper 
Lx — Tx relation. Observations suggest that such steepen ing 
does happen below 1 keV fe.g.. iHelsdon fc PonmanlboOOl ). 

More quantitatively, we directly calculate the Lx — Tx 
relation using the density and temperature profiles from our 
one-dimensional models (e.g., see Fig. [T}. We calculate the 
luminosity-weighted temperature 

_ n £ n lT A(T)4nr 2 dr 

X j; 200 n e n l A(T)4 7 rr 2 dr ' 1 ' 

and X-ray luminosity 

/•1-200 

Lx = I n e rii A{T)4-Krdr. (4) 
Jo 

The solid line in Figure [3] shows the Lx — Tx relation for our 
fiducial model with an entropy core and tTi/tg ^ 10 (corre- 
sponding to solid lines in Fig. [1]); the asymptotic gas density 
slope is —2.25. The Lx — Tx relation for the fiducial case pro- 
vides a very good fit to the upper envelope of observational 
points. This is very satisfying because our tn/iff = 10 model 
corresponds to an upper limit on the core density and hence 
on Lx ; the core density and luminosity is lower than this for 
non-cool core halos. The models with a constant tTi/iff in the 
core, both with and without the BCG potential, give results 
similar to the entropy core models because the contribution 
of the increased density in the core is small. 

The long-dashed line and the dot-dashed line in Figure 
[3] show the core entropy models with tTi/tg = 40, 100 re- 
spectively. For massive halos, the Lx — Tx relation does not 
depend sensitively on the critical tTi/ts used for introducing 
the entropy core because the cores are small in these halos 
and contribute a subdominant fraction of the total luminos- 
ity. In lower mass halos, the core sizes are larger; thus, Lx is 
more sensitive to the tTi/tg threshold. Any 'overheating' by 
AGN or supernova feedback (which is reflected as an increase 
in tTi/tg) could therefore contribute to the spread observed 
in the Lx — Tx relation for lower mass halos. 

Figure [3] also shows that the outer gas density slope 
(s = dlnp/dhir) affects the Lx — Tx relation for the high and 
low mass halos differently. The luminosity is higher (lower) 
for clusters (groups) with steeper outer density profiles than 
for clusters (groups) with shallower profiles. For our lowest 
mass halos, where tTi/ts < 10 even at r2oo, the halo pro- 
files are expected to be similar. The X-ray luminosity for a 
density profile with dlnp/dlnr — —3 is ~ 4 times larger 
than for dlnp/dlnr = —1.75 at the massive halo end. Thus, 
the non-universality of the outer gas density profile (e.g., 
ICroston et al. 1 12008). along with the variation in the core 
density, can explain the large spread in Lx — Tx relation for 
massive halos. 

Recently iMaughan et al.~l (|2012T ) have analyzed clusters 
with different temperatures and concluded that the relaxed 
hot (j> 3.5 keV) clusters have self-similar densities. Corre- 
spondingly, they find these clusters to satisfy the self-similar 
Lx — Tx relation, Lx oc T x - The solid line in Figure [3] shows 
that this is also the case for Tx J> 1 keV clusters in our fidu- 
cial model corresponding to cool core halos. Self-similarity 
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3. The X-ray luminosity (Lx', Eq. SJ as a func- 
the luminosity- weighted temperature (Tx', Eq. |3j 
using our models and fr om observations. The observ a tiona l 
data points are from [jeltema, Binder, fc Mulchaey I [200; 
Holden et al. Il200l iMarkevitch Ill998l. iHelsdc-n fc Ponman ibopl 
"20031, iDella Ceca et al. Il2000l, lArnaud fc Evrardl 
2009, and references 
therein. The relations obtained from our one-dimensional models 
are: entropy core models with frriAff 10 and different outer 
gas density slopes (s = dlnn/dlnr); entropy core models with 
tTl/tff ^ 40, 100 and the outer density slope of —2.25. Increasing 
the tii/iff threshold reduces the X-ray luminosity, especially for 
low mass halos. 



Mulchaey et al. .._ 
19991 . iBorgani et al. I l200ll. IPratt et al 



seems to persist until 1 keV in Figure [3] but in reality the 
density in smaller halos may be lower than given by Figure 
Q] because of over-effectiveness of feedback in low- mass halos 
(Fig.|4]shows that the required feedback is smaller than AGN 
feedback efficiency for halos masses <C fewxlO 14 Mg). 

An important fact that must be noted when compar- 
ing our models to observations is that the data in Figure [3] 
are from different cluster/group/galaxy samples, both local 
and cosmological, and using different X-ray instruments. A 
lot of these observations are biased. For example, low tem- 
perature groups and galaxies have low surface brightness at 
large radii such that significant X-ray lum i nosity from close 
to the virial radius is missed (e.g., |V02| . iMulchaev et al. 1 
l200l iHelsdon fc Ponman I l200d i. Some of the data points 
have the contrib ution from the cool core removed (e.g., 
lMarkevitch~lll99ct ). and some samples only con tain clusters 
witho ut cool cores (e.g., Arnaud fc Evrardll 19991 ) . The galax- 
ies in ljeltema. Binder, fc Mulchaey 1 (|2008h are from groups 
and clusters, and are probably fainter than the field galax- 
ies because of ram pressure and tidal stripping. Some data 
are from clusters at cosmological distances, thus its more 



difficult to detect plasma close to the vi rial radius (e.g., 
IDella Ceca et al. II2000J ; [Holden et al. 1120021 ). 

While we calculate Lx out to the virial radius T2oo , most 
observations only go out to rsoo- Recent observations show 
that the baryon content becomes closer to the universal value 
if the baryons are coun ted out to 7-200 (e.g.. lSimionescu et al. I 
l201ll ; ISato et al.~ll2012l ). Thus, for lower mass halos, where the 
density profile is shallower (see Fig. [TJ and surface bright- 
ness is lower, large radii contribute significantly to Lx, and 
the observed luminosity is expected to be smaller than pre- 
dictions. All these uncertainties preclude us from getting an 
exact match with observations. The overall trend of observa- 
tions compared wi th predictions is very satisfying, however. 
ICrain et al. (2010) argue that the Lx — Tx relation for the 
X-ray halos of spiral galaxies is similar to the relation for el- 
lipticals. This is expected if the X-ray halo is maintained in 
tcooi/tff > 10 state in both spiral and elliptical galaxies via 
supernova and AGN feedback, respectively. 

Previous works have tried to explain the low luminosity 
of galaxy groups relative to the scaling relations by expul- 
sion of gas from smal l er halos through AGN feedback (e.g., 



Puchwein et al. 2008); iMcCarthv et al. Il20ld : iFabian et ail 
2010f ) or by mass dropout from the hot ph ase into the cold 



phase (e.g.. lBrvanll2000l : iNaeal et al. 1 I2007I ). However, there 
are problems with both scenarios. The mass dropout scenario 
without sufficient heating predicts massive cooling flows and 
excessive star formation, not observed even at group scales. 
The scenario where a large fraction of the X-ray emitting gas 
is blown out beyond the virial radius for groups, but not for 
clusters, is at odds with observations which show that the 
total halo baryon fraction within rsoo, including the contri- 
butio n of intracluster stars , is similar for clusters and groups 
(e.g., iGonzalez et al. II2007T I. The gas fraction decreases with 
a decreasing halo mass, but the decrease in gas fraction is 
roughly compensated by the increase in stellar fraction. By 
modeling stellar fraction in more detail, rec ent papers, e.g., 
iGiodini et al~ll2009l ; iLeauthaud et al. II2OI2L show that there 
is a slight (less than a factor of 2) decline in the baryon frac- 
tion on going from massive groups to 10 13 Mq halos. However, 
even for 10 13 Mq halos the ICM mass is larger than stars, and 
a smaller X-ray surface brightness prevents an accurate mea- 
surement of the gas fraction, especially out to r2oo- 

Our models roughly reproduce the observed den- 
sity/entropy profiles (Fig. [1} and the observed Lx — Tx rela- 
tion (Fig. [3}. Moreover, since cooling (and probably feedback; 
sec Section |3~2]FI is not important close to the virial radius for 
halos > 10 13 Mq (Fig. [TJ, we expect groups and clusters to 
be approximately closed systems (see fd in Fig. [4] and Table 
HJ, as the observations of IGonzalez et al.~l (|2007l ) suggest. 



2 Even if cooling is not important close to the virial radius, heat- 
ing/outflows can, in principle, transport baryons out of the virial 
radiu s and thus reduce the baryon fraction (e.g.. [Puchwein et al. I 
120081 ). 
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Figure 4. Various quantities as a function of the halo mass for 
our fiducial model ( results are not too different for in Aft = 10 
'core' model): the efficiency required to balance cooling (e req = 
Lx /MrfC 2 ; see Table[l]for the definition of various quantities); the 
mass fraction expected to be present in the ICM, /icm ~ (1 — fd)fb\ 
the dropout mass fr action, fdfb \ fit t o the baryon mass fraction, 
. using Eq. 7 from iDai et al. | |201C ); and the stellar mass frac- 
tion, /*, from Figure 5 of lDaiet^h* |2010| . Unlike them, we use 
the halo mass M200 as the independent variable, instead of the cir- 
cular velocity, using Af 2 00 = 10 15 M [V c /1632.07 km s -1 ] 3 . Also 
indicated are the typical ICM heating efficiencies expected due to 
AGN and supernovae (ebh an d e*). Notice the sudden decrease 
in the required efficiency and the change in other quantities at 
~ 10 13 Mq. For halos < 10 13 Mq we expect most baryons to be 
beyond r200 and only a small fraction (/„ given by Eq. \E\ green 
long-short-dashed line) in form of stars. 



3.2 Baryon Budget in Halos 

In this Section we use our models to study the baryon frac- 
tion as a function of halo mass. In our models the cluster 
mass halos tend to develop small cores (Fig. [T]) and thus the 
dropout mass is insignificant; thus, clusters are baryonically 
closed with a baryon fraction similar to the cosmic value. 
Lower mass systems, however, can have a lower baryon frac- 
tion. In fact, for halo mass < 10 13 Mq, iTi/iff < 10 even at 
the outer boundary (see Fig. QJ and the ICM is depleted of 
most of the baryons, either due to star formation or because 
of expulsion of gas beyond the virial radius by feedback heat- 
ing (see f d in Table [T]). Our models do not directly predict the 
fate of baryons absent from the ICM: whether they form stars 
or are lost due to feedback. However, the efficiency required to 
balance cooling (e rcq = Lx/Mjc 2 ; see Table [TJ as a function 
of halo mass provides some insight. If the required efficiency 
is small compared to the available feedback source, we expect 
the halo gas to be overheated, and sometimes even pushed be- 



yond the virial radius. This is the case for halos less massive 
than 10 13 Mq. Thus, we expect these halos to have substan- 
tial baryonic mass 'missing' from within r2oo- This feedback- 
heated hot gas beyond the virial radius is different from the 
warm-hot intergalactic medium (WHIM; iDave et al. II2OO1T ) 
which is heated by structure formation shocks in cosmologi- 
cal filaments. 

As argued above, a closed-box model (in which the inflow 
has the universal baryon fraction and outflows are negligible) 
for baryons is not good for halos less massive than 10 13 Mq. 
The equations of mass and energy conservation in halos, in 
steady state, are 



Min ~ Micm + Mbh + M* + Mo 



(5) 



e*M*c 2 + 6bhMbhc 2 w Lx.iCM+out + Excm + E ont , (6) 

where Mi n is the baryonic mass accretion rate onto the 
halo, AficM is the rate of increase in the ICM mass, Mbh 
is the mass accretion rate onto the central SMBH, M* is 
the star formation rate, and M out is the mass outflow rate 
(outflows are driven by supernova and AGN feedback). Sim- 
ilarly, e„(eBH) is the supernova (AGN) efficiency in heat- 
ing the ICM, Lx,iCM+out is the X-ray luminosity of the 
ICM and the hot gas beyond T2oo, -Ercijfl is the rate of in- 
crease of the total ICM energy (internal + gravitational), 
and E out w (5/2) f v irM out k B T2oo/ nm v (f vll is a factor com- 
paring the temperature of the outflowing gas to X200) is 
the total energy outflow rate from the halo. Eqs. 1 5161 are 
severely underconstrained as stated, but observational data 
and our models can be used to constrain parameters such 
as Af„/Min, Mbh /M in , M out /M in . For example, the dropout 
fraction fd = Md/M[ n ~ 1 — Micm/ Mi n can be read off from 
Table [T] In massive halos (> 10 13 M and f d <C 1), the dom- 
inant balance in mass conservation is Mi n « Micm; i.e., most 
of the mass accreted by the halos is incorporated in the ICM. 
Similarly, for massive halos Eq. [6]reduces to ebhMbhc 2 ~ Lx 
because supernova feedback is expected to be subdominant, 
and £icm and _E out are negligible □ 

Figure|4]shows the feedback efficiency required to balance 
cooling, e rcq = Lx /Md(? , as a function of the halo mass. The 
required efficiency decreases w it h a decreasing halo mass, as 
noted by iGaspari et al. I (|201ll ) ; iGaspari, Brighenti, fc Temil 
(|2012l ). The halos can be divided into three categories, de- 
pending on e r cq: first, the most massive halos (> 6xl0 14 M Q ) 



3 Note that in absence of feedback heating and cooling, Eq. [6] 
becomes -Bicm ~ 0, which reflects rough virial equilibrium. The 
weakly bound inflow virializes at the accretion shock and the grav- 
itational energy is converted into thermal energy. Moreover, the 
timescale corresponding to -Eicm is the halo age, which is much 
longer than the cooling time for groups and clusters. Detailed mod- 
eling of -Bicm is beyond the scope of this paper. 

4 Note that a substantial contribution to Lx comes from outer 
radii where the cooling time is long. Thus, AGN feedback need not 
power all of Lx', some of Lx can come at the expense of -Eicm- 
Similarly, not all of the dropout mass is accreted by the SMBH. 
Here, we simply equate Lx and CreqAf^c 2 to avoid introducing 
more free parameters. 



© 2012 RAS, MNRAS 000.flM3l 



Lx — Tx Relation & Missing Baryons 9 



in which e rcq is larger than the efficiency of kinetic /radio- 
mode AGN feedback (e B H ~ 0.003; e.g., iMerloni fc Heinz I 
120081 ); second, the intermediate mass halos (10 13 - 10 14 M Q ) 
in which e roq lies between 6bh and e» (~ 10" 6 - 1O" 7 )0 and 
third, the lowest mass halos in our study 10 13 M Q ) m 
which e rcq <C e» . 

After introducing the three halo classes based on e rcq , 
we can study heating in each of these classes. The required 
efficiency for the most massive halos (~ 10 15 M ) is > 0.003 
(see Table [1]), larger than typical estimates for radio mode 
AGN feedback (~ 0.003), especially given the fact that only 
a small fraction of dropout mass may feed the SMBH. Thus, 
a substantial contribution to heating in massive halos is prob- 
ably due t o non- feedback processes such as thermal conduc- 
tion (e.g., iGuo. Oh. fc Ru szkowski~1 l2008h . For intermediate 
mass halos (10 13 - 10 15 M Q ) AGN feedback is adequate to 
balance cooling (Fig. [4j; however, supernova feedback is still 
insufficient. While the required feedback efficiency (e rcq ) is 
smaller than ebh, the majority of dropout mass may lead to 
star formation instead of accretion. However, if a fixed frac- 
tion of Md is channeled into Mbh for all halos, galaxy groups 
are expected to be mostly in a non-cool core state because 
a larger feedback efficiency res ults in a lower frequency of 
cooling episodes (see Fig. 11 in iPaper it ). For intermediate 
(10 13 — 10 15 Mq) halos feedback heating is unlikely to ex- 
pel baryons beyond the virial radius because the core is well 
within the virial radius (Fig. [TJ and excess heating is ex- 
pended in expansion of the core (Eicm in Eq. [6j . Thus, even 
intermediate class halos are expected to be roughly baryoni- 
cally closed. 

The required efficiency (e rcq ) is very small (< 10~ 7 ) for 
the lowest mass halos (<; 10 13 M ; see Fig. [IJ. The feed- 
back efficiency expected from star formation (mainly via su- 
pernovae) is ~ 10 -6 — 10 -7 , smaller compared to the AGN 
feedback efficiency. But even this is much larger than the effi- 
ciency required to power the X-ray luminosity in lower mass 
halos (see Fig. 3]). More quantitatively, a comparison of e roq 
and e» in Figure [4] shows that every baryon in the hot phase 
of smaller halos (< 10 13 Mq) has ~ 10 — 100 times more en- 
ergy available than is required to power its X-ray luminosity. 
Because of such a strong feedback heating, we expect most 
dropout baryons in lower mass halos to not form stars but 
instead be expelled beyond the virial radius. Moreover, the 
fraction of 'missing' baryons in low mass halos is expected to 
increase with the decreasing halo mass because of a decreasing 

^rc q • 

Assume that the fraction of mass accreting onto the 

5 We estimate the efficiency of supernova feedback as follows. We 
assume that th e stars with 8 Mp) < M« < 20 M become su- 
pernovac (e.g., IWooslev. Heger. fc Weaver 1120021) and assume a 
Salpeter IMF. The mass fraction in stars which will become super- 
novae is fa 0.07. Only a small fraction (~ 0.1) of the total supernova 
kinetic energy (~ 10 51 erg) couples to th e hot phase because of ra- 
diative losses (e.g.. lThornton et al. Il998h . Thus, the supernova effi- 
ciency for heating the ICM is ~ O.O7xO.lxlO 51 /8M c 2 ~ 5xl0 -7 
(indicated in Fig. [4}. 



SMBH in low mass halos is negligible. Let us define /« = 
f b M*/M in , /out = /bMout/Min, /icm = j&MicM/Min, and 
e req = Lx/MdC 2 , where M d = M* + M out . Thus, from Eq. 
El fbfd ~ /* + /out and /icm + /out + /» « fb- The energy 
conservation equation (Eq.[6j becomes 

uM*c « L x 1 + -, + /virMout, (7) 

V /icm / 2 fj,m p 

where supernova feedback is assumed to power the X-ray lu- 
minosity of the ICM and the extended halo, and the thermal 
heating of the outflowing hot gas. Using our definitions, the 
above equation reduces to 

f* w 1 e. + (5/2)/ vir (l - f d )k B T 2Q0 /fim p c 2 

fb ~ [e» + £req/d/(l - fd) + (5/2) f vil k B T 2 oo / fJ,m p c 2 ] ' 

(8) 

where / v i r is a virial factor comparing the temperature of the 
outflowing gas to the virial temperature (T^oo)- Note that the 
above equation is only valid when e* > e roq , i.e., only for halos 
<C 6 x 10 Mq . Eq. [8] should only be treated as an order of 
magnitude expression because the cooling time for low mass 
halos at large radii is long (~ 10 Gyr) and our assumption of 
thermal balance need not hold as tightly. Moreover, the frac- 
tion of the supernova feedback energy going into heating the 
outflowing gas (the / v i r terms in Eqs. [7] & [HJ is very uncer- 
tain. The exact distribution of baryons in stars vs. 'missing' 
baryons (beyond the virial radius) will require detailed nu- 
merical simulations beyond the scope of this paper. Here we 
only provide an order-of-magnitude estimate for the baryon 
budget in low mass halos. 

Figure U shows /» as a function of the halo mass for lower 
mass halos using Eq. [8j and setting / v i r to zero (green long- 
short-dashed line). The expectation is that the overestimated 
X-ray luminosity in Eq. [7] (its an overestimate because the 
density beyond the virial radius is assumed to be the same as 
the core density) is of the order of the contribution of feedback 
in driving outflows. Plugging in / v i r < 1 in Eq. [HJ should not 
change our results significantly. The predicted stellar frac- 
tion for low mass halos shows a good match with observa- 
tions. Figure U also indicates that excessive feedback in halos 
< 10 13 Mq results in thermal expulsion of majority of baryons 
beyond the virial radius (/ out w fb — /* - /icm > /*, /icm). 
It is likely that these baryons were never incorporated within 
halos but were kept beyond r2oo bec ause of excessive heat- 
ing, a s suggested observationally by lAnderson fc Bregman I 
([20101 ). 

Recall that we choose ft = 0.17, the universal value, in 
constructing our models (see Table [TJ. However, as we argued 
in the previous paragraph, most baryons in low mass halos 
(< 10 13 Mq) are likely to be expelled beyond r2oo- Thus, 
the baryon fraction for lower mass halos, only accounting for 
baryons within the virial radius, should be much smaller than 
the universal value. Indeed, the total baryon fraction within 
f"2uo, /* + /icm, in Figure [4] is much smaller than 0.17 for 
halos less massive than 10 13 Mq. Table[JJshows that the mass 
dropout rate (M d ) for the Milky Way halo (10 12 M ) is w 
28.2 M yr~\ but our models suggest that the majority (~ 
90%) of these baryons are likely to be beyond r^oo (Fig. HJ . 
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Thus, the star formation rate expected for the Galaxy is ~ 
O.lMd, which is ~ 2 — 3 Mp yr -1 (see Tabled]), quite sim ilar 
to the observed value l|Bauermeister. Blitz, fc Mall2O10l and 
references therein). 

Observationally, most of the ba ryons for halos > 10 1 3 Mq 
are in the form of the ICM (see also iGiodini et al. H2009I 1 and 
a closed box model for massive halos is a good approxima- 
tion. Our prediction for the mass fraction in the ICM (dotted 
line in Fig . 2| agrees qualitatively with the best-fit line from 
iDai et ahl (|2010l ). In particular, the sharp decline in baryon 
fraction for halos <J 10 13 Mq ( corresponding to the temper- 
ature of ~ 1 keV, as noted bv lDai et al. II201QT ) is correctly 
captured by our models. Quantitative agreement is difficult 
because observational data show large scatter about the best 
fit. The observed mass fraction in stars also shows a large 
scatter and is much larger than the predicted dropout frac- 
tion for massive halos 10 13 M ) because our models only 
consider accretion of baryons in form of gas. However, galax- 
ies composed of stars, falling along cosmological filaments, are 
also accreted by massive halos. Thus, Figure [4] suggests that 
the stars from previous generation of galaxies (which merge 
to form the central BCG, and contribute to other galaxies 
and intracluster stars) form the majority of stars in galaxy 
clusters; the contribution of stars due to cooling of the hot 
ICM is negligible. 

Galaxy clusters seem to pose an apparent paradox. The 
observed stellar mass fraction in clusters is smaller than in 
10 12 — 10 13 Mq halos (see Fig. 3}. However, clusters are formed 
via hierarchical mergers of star-rich smaller halos. How can 
the stellar fraction in clusters be smaller than in smaller ha- 
los? This paradox can be resolved by noting that the stellar 
fraction in even smaller (-C 10 12 Mq) halos falls precip itously 
with a decreasing halo mass (Fig. 5 in lDai et all2010f ). Thus, 
if a significant fractio n of cluster growth occurs via accretion 
of very small halos (jFakhouri. Ma. fc Bovlan-Kolchinl l20ld 
suggest that smooth accretion contributes predominantly to 
the halo growth), clusters can have a smaller stellar fraction 
than ~ 10 12 Mq halos. As discussed in the previous para- 
graph, our models which only consider gas accretion obtain 
a stellar fraction much lower than in observed clusters. The 
detailed prediction of stellar fraction for clusters will require 
halo merger trees, and is beyond the scope of the present 
paper. 



4 CONCLUSIONS 

Inspired by our previous numerical simulations of the local 
thermal instability in the ICM with feedback heating, we have 
constructed simple one-dimensional models of the ICM. Ac- 
cording to our simulations, local thermal instability in the 
ICM leads to multiphase gas if, and only if, the ratio of the 
thermal instability timescale and the free-fall time (in /iff) is 
< 10. For densities higher than this (i.e., t-ri/iff < 10) mas- 
sive cold filaments will condense out of the hot phase, and the 
hot ICM will adjust roughly to a state with in/iff ^> 10 ev- 
erywhere. Figure Q] shows that the core density decreases with 



a decreasing halo mass; the decrease is quite pronounced for 
halos less massive than 10 13 Mq. Observations also show that 
the c ores are bi gger and more tenuous in smaller mass halos 
(e.g.. lSun~ll2012r i. A lower core density for smaller halos natu- 
rally produces the steepening of the Lx —Tx relation for lower 
mass halos (see Fig- EJ) - We do not expect the density profile, 
the Lx — Tx relation, and the baryon budget to depend sen- 
sitively on redshift, as long as the feedback sources (AGN, 
supernovae) are active and the cooling time is shorter than 
the halo age. Thus, the assumption of a constant gas fraction 
for clusters at all redshi fts, which is used t o constrain cosmo- 
logical parameters (e.g.. lAllen et al."ll2008r ). is expected to be 
valid. 

According to our models, the density of the Galactic halo 
is expected to be small (n e ~ 10~ 4 cm -3 ) because of a short 
cooling time at the virial temperature of the Galactic halo 
(~ few x 10 6 K). Thus most of the plasma in the ICM is 
expected to either cool and form stars or be expelled beyond 
the virial radius (see fd in TableQJ. Due to excessive feedback 
heating by supernovae, most of the baryons are pushed out of 
the virial radius instead of forming stars (see /» and /icm for 
halos < 10 13 Mq in Fig. [J). Our models show that the halos 
less massive than 10 13 M should be 'missing' a substantial 
fraction of their baryons, in agreement with observations (e.g., 
IDai et al. Ikoiol ). 

We tested two different models for the central cores, 
which make very different predictions for the gas density 
in the Galactic halo (~ 10 12 Mq). Observationally it has 
not been possible to unambiguously characterize the Galac- 
tic h alo because of the low number d ensity in the hot phase 
(e.g.. [Bregman fc Llovd-Davis I [20071 ). However, pulsar dis- 
persion measure and absorption line studies suggest that 
the number density is quit e low (n e ~ 10 -4 cm -3 ; e.g., 
lAnderson fc Bregman II2OI0I ) . A recent study which combines 
the Galactic absorption and emi ssion observations of the hot 
gas points to a similar value dGupta et al. I l2012f ). There 
are also indications (e.g.. lEverett et al. 1120081 and references 
therein) that some portions of the Galactic hot gas have den- 
sities as high as ~ 0.01 cm -3 , consistent with the prediction 
of our constant txi/tg 'core' models which include the grav- 
itational influence of the central galaxy. 

Extended X-ray emission correlated with the indicators 
of star formation (UV, FIR, Ha) has been observed in ha- 
los of some spiral galaxies. The X-ra y emission is believed 
to be powered by su pernovae (e.g., iTiillman et al. 1 120061 ; 
IStrickland et al. 1120041) but some of it may also be due to 
the larger emissivity of a denser ICM. If so, this correlation 
of the low-entropy, high emissivity X-ray gas and the sig- 
natures of feedback (star formation in case of halos of spi- 
ral galaxies) is analogous to the correlation between the en- 
tropy and cold gas/AG N feedback in cluster s of galaxies (e.g., 
ICavagnolo et al.~ll2008l ). Along similar lines. iTumlinson et al. I 
( 20111 ) have recently observed large quantities of cooling gas 
(OVI) in outskirts of star forming galaxies, indicating a causal 
connection between the cooling of the hot halo and star for- 
mation. Our thermal instability + feedback models for hot 
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halos try to explain the baryonic structure of all halos — from 
galaxies to clusters — in a single framework. 

According to our nonlinear criterion for the condensation 
of cold filaments in a gravitational field (frn/iff < 10) the high 
velocity clouds (HVCs) of atomic hydrogen in Galactic mass 
halos can condense out of the hot phase close to (or even 
beyond) the virial radius (~ 200 kpc). This multiphase gas 
at distances large compared to the optical size of galaxies 
can explain quasar absor ption lines observ ed in the neigh- 
borhood of galaxies fe.g.. IChen et ain i2010). A much larger 
density of the hot gas is required for the HVCs to condense 
out at smaller radii because of a short free-fall time (see the 
dotted line corresponding to the iTi/iff = 10 'core' model 
for the halo mass of 10 12 M in Fig. [TJ|. Once HVCs form, 
these dense, low metallicity filaments will fall freely toward 
th e center of the halo, as observed (e.g., observations of M31 
by iThilker et al~ll2004l ) . These clouds will not only be shred- 
ded by Kelvin-Helm holz instabilities at the interface (e.g., 
iMurrav et al~lll993l ) but also reform due to efficient cooling 
of the shredded cool gas (see Fig. 4 in lPaper III) . 

We have intentionally kept our models simple by using 
the NFW gravitational potential for most halos; we have 
also included a fixed BCG potential for some of the mod- 
els (see Fig. [2]). The BCG models can explain the non-cool 
core clusters and groups which show much small er but dens e 
'coronae' that can power AGN feedback (e.g., ISunll2009h . 
Similarly, some isolated elliptical galaxies, in which the stel- 
lar bulge contributes significantly to the gravitational poten- 
tial at inner radii and reduce s the free-fall time, can have 
larger central densities (e.g ., IWerner. Allen, fc Simionescu I 
120121 ; iHumphrev et al. II2012I ) . 

From our models based on tTi/tff, we expect the ellip- 
tical galaxies in cluster centers (the BCG, in particular) to 
have gas densities much larger than similar mass field galaxies 
because the virial temperature, and hence the cooling time, 
is smaller for isolated galaxies. There is some observational 
support for this hypothesis (e.g.. iMulchaev fe Jeltema~ll2010l ; 
ISun et al. II2009I ). but more observations are needed. We ex- 
pect a lower ICM density for non-central elliptical galaxies in 
clusters, where the gravity of the dark matter halo does not 
point toward the center of the elliptical galaxy; ram pressure 
and tidal stripping introduce additional effects which further 
lower the hot gas density. Detailed environmental effects on 
the structure of hot gas in elliptical galaxies can only be stud- 
ied with numerical simulations, and not by our simple models. 
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Table 1. One Dimensional Models with an Entropy (txiAff) Core using critical iTiAff = 10 and s = —2.25 



Mjoo ( M o) 


ct 


n e ,o (cm 3 ) 


Kq (keVcm 2 ) 


T x (keV) 


L x (xlO^ergs" 1 ) 


ff 


M'A M yr" 1 ) 


e tt q 


3 x 10 15 


4 


0.098 (0.15) 


20.7 (11.1) 


13.8 (13.7) 


79 (79.2) 


2.88 (0.37) x 10 -4 


29.3 (3.8) 


0.0047 (0.037) 


10 15 


5 


0.071 (0.13) 


22.6 (8.68) 


7.07 (7.04) 


24.1 (24.2) 


4.14 (1.1) x 10~ 4 


14.1 (3.75) 


0.003 (0.013) 


6 x 10 14 


6 


0.069 (0.14) 


22.3 (7.7) 


5.4 (5.3) 


15 (15.2) 


4.82 (1.84) x l(r 4 


9.8 (3.75) 


0.0027 (0.0071) 


3 x 10" 


7 


0.058 (0.13) 


22 (6.58) 


3.6 (3.59) 


7.3 (7.44) 


6.12 (3.8) x 10~ 4 


6.25 (3.9) 


0.0021 (0.0034) 


10 14 


8 


0.038 (0.095) 


21.3 (5.3) 


1.91 (1.88) 


2.1 (2.17) 


1.32 (1.31) x l(r 3 


4.5 (4.4) 


0.00082 (0.00085) 


6 x 10 13 


9 


0.032 (0.088) 


21 (4.77) 


1.46 (1.44) 


1.23 (1.28) 


2.77 (2.37) x l(r 3 


5.65 (4.84) 


0.00039 (0.00046) 


3 x 10 13 


9.5 


0.02 (0.065) 


23.2 (4.34) 


0.98 (0.96) 


0.51 (0.55) 


8.8 (5.8) x 10" 3 


8.99 (5.95) 


9.86 x 10" 5 (0.00016) 


10 13 


10.5 


0.0045 (0.035) 


47.9 (4.04) 


0.5 (0.51) 


0.078 (0.13) 


0.08 (0.036) 


27.3 (12.1) 


5 x lO -6 (1.88 x 10" 5 ) 


6 x 10 12 


11 


6.4 x 10" 4 (0.023) 


159.6 (4.11) 


0.36 (0.38) 


6.79 x 10" 3 (0.045) 


0.5 (0.2) 


102.6 (41.3) 


1.15 x 10 -7 (1.92 x 10 -8 ) 


3 x 10 12 


11.5 


3.3 x 10~ 4 (0.011) 


171.1 (4.56) 


0.27 (0.26) 


1.58 (7.7) x 10" 3 


0.68 (0.52) 


68.9 (53.4) 


4 x 10" 8 (2.5 x 10~ 7 ) 


10 12 


12 


1.2 x 10" 4 (3 x 10" 3 ) 


193.3 (6.16) 


0.18 (0.16) 


1.88 (5.69) x 10- 4 


0.83 (0.78) 


28.2 (26.4) 


1.2 x llr 8 (3.8 x 10~ 8 ) 



The halo baryon mass fraction is chosen to be the universal baryon fraction /j, = 0.17 for all halos. The entropy core model is the fiducial 
model. The quantities in parentheses denote the values for a constant txiAff 'core' model. The various quantities are: n e ,o the central 
electron density (at r = 0.001r2oo)> the central entropy, Tx the luminosity-weighted temperature (Eq. [3j , and Lx the X-ray luminosity 
(Eq.H). 

^Af200 is the halo mass within r200i such that, M200 = (47r/3)200p C r^200' w h ere Per is the critical density taken to be 9.2 X 10 — 30 g cm -3 . 

ic is the concentration parameter defined as c = r^ooAs, where r s is the scale radius. 

tt/^ = (/i,M2oo — Jc^ 200 ^ 7rr2 pdr) I fbM200 is the mass dropout fraction, where p is the gas mass density. 

*Md = fbfd Af200 /*age ~ M ou t + M, + AfflH ~ Afi n — A/tcm ( see Eq. O is the mass dropout rate, where we choose t age = 5 Gyr. 
TTe req = Lx/Md c2 is the heating efficiency required to balance radiative cooling of the ICM. 
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